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ABSTRACT 

We present medium resolution optical spectroscopy of a sample of nine Low Surface 
Brightness (LSB) galaxies. For those that show clear signatures of AGN emission, 
we have disentangled the AGN component from stellar light and any Fe i and Fe ii 
contribution. We have decomposed the Ha line into narrow and broad components 
and determined the velocities of the broad components; typical values lie between 
900-2500 kms '. Of the galaxies in our study, UGC 6614, UGC 1922, UGC 6968 
and LSBC F568-6 (Malin 2) show clear signatures of AGN activity. We have calcu- 
lated the approximate black hole masses for these galaxies from the Ha Une emission 
using the virial approximation. The black hole masses are ~ 3 x 10^ Mq for three 
galaxies and lie in the intermediate mass black holes domain rather than the super- 
massive range. UGC 6614 harbors a BH of mass 3.8 x 10^ M©; it also shows an 
interesting feature blueward of Ha and H/3 implying outflow of gas or a one-sided jet 
streaming towards us. We have also measured the bulge stellar velocity dispersions 
using the Ca ii Triplet lines and plotted these galaxies on the M - cr plot. We find 
that aU the three galaxies UGC 6614, UGC 6968 and F568-6 lie below the M - cr 
relation for nearby galaxies. Thus we find that although the bulges of LSB galaxies 
may be well evolved, their nuclear black hole masses are lower than those found in 
bright galaxies and lie off'set from the M - a correlation. 
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1 INTRODUCTION 



Low S urface Brightness (LSB) galaxies are an extreme class of late type spiral ga laxies (|Impey & Bothun 



19971). 



hey are poor in star formation (IBoissier et al 



licities (IMcGaugh 



2008 



O'Neil et al 



1994 ) and low dust masses (IRahman et al. 



20071 : 



2004h, have low metal- 



Hinz et al. 



2MB- However, 



they are gas rich (H i) and their Hi disks are far more extended than their diffuse stellar disks 



( O'Neil etal. 



20041). Their disks have weak spiral arms and their bar perturbations are not as 



prominent as those seen in bright galaxies. This low level of disk activity can be attributed to 



the presence of dark matter ha' 



instabilities (iMayer & Wadsley 



OS (de Blok et al. 



2004; 



Mihos et al. 



2001 ) that tend to prevent the formation of disk 



1997 



Galaz et al. 



201 1|) . Surveys show that al- 



though LSB g alaxies are common in ou r local universe they are preferentially found in low density 



environments (Rosenbaum et al. 



20091) and can have wide ly varying morphologies; starting from 



the populous dwarf LSB galaxies (ISchombert et 



which Malin 1 is a good example (|Pickering et al 



20011) to the giant LSB galaxies (GLSB), of 



1997). 



Although a lot of work has been done towards understanding the stellar and gas content of 
LSB galaxies, not much is known about their nuclear activity or central black hole (BH) masses. 
A significant fraction of G LSB galaxies hav e active galactic nu clei (AGN) that are often asso- 



ciated with massi ve bulges (|Schombert 



(Das et al. 



1998 



Im pey et al.ll200lh. The AGN may be radio bri ght 



2009ah and visible at X-ray wavelengths as well (|Naik et al. 



201 



Das et al. 



2009b). It 



is now well established that the black hole masses in gal axies are correlated with their bulge veloc 



ity di spersions (M - cr) and bulge luminosities (M - L) (|Ferrarese & Merritt 



20001 : 



Gebhardt et al. 



2000() wh ich suggests that black hole formation, galaxy evolution and AGN act ivity are all in 



20101 : 



Somerville et al. 



2008 



Merloni & Heinz 



20081). However, the 



terlinked (|Schawinski et al.l 
M - cr a nd M - L relations s how less scatter when only ellipticals and early type galaxies are 



included (Gtiltekin et al. 



2009h . Late type spirals introduce more scatter in the correlations. This 



may be due to a larger intrinsic scatter in the BH masses of late type spirals or simply measurement 
uncertainities since late type spirals, in genera l, have srn aller bulges. They also show weaker AGN 



activit y compared to the earlier type galaxies 



Beifiori et al. 



HI2QO8D. Alternatively, In a study by 
(120091) involving HST observations of 105 nearby galaxies spanning a wide range of Hubble types 
from ellipticals to late-type spirals, the estimated black hole masses upper limits appear to lie 



closer to the expected black hole masses in the most massive elliptical galaxies with values of cr 
above 220 km/s than for galaxies with cr in the range 90-220 km s"' which appears to be consistent 
with a coevolution of supermassive black holes and galaxies driven by dry mergers. 
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It is not clear where LSB galaxies lie on the M -a plot. Their bulge velo city and disk rotatio n 



speeds suggest that they lie below the M - cr correlation for bright galaxies (jPizzella et al 



2005) 



X-ray studies also suggest that GLSB galaxies do not lie on the radio-x -ray correlation (|Das et al 



2009b|) and their black hole masses maybe quite low (|Naik et al.ll2010l) . In this paper we use de 



tailed optical spectroscopy to determine the position of a sample of LSB galaxies on the M - cr 
diagram. We observed the nuclear spectra of several bulge dominated LSB galaxies; for galaxies 
that showed AGN emission, we estimated both the black hole masses and also the bulge veloc- 
ity dispersion. In the following sections we present our observations, the results and discuss the 
implications of our findings. 



2 SAMPLE SELECTION 



Our s ample consists of nine large, bulge dominated 



19871), of which eight have been observed further by 



.SB galaxies from 



Schombert & Bothun 



Schomberd (|l998|) . The basic parameters 



of the galaxies are listed in Table [B The LSB galaxies in the Schombert sample were all Hi rich, 
giant, spiral galaxies that were derived from the UGC catalogue; they have systemic velocities that 
are less than 15, 000 km s"' . We chose a subset of eight nearby galaxies from that sample that had 
Vsys < 10, 000 km s"' and appear to be bulge dominate d galaxies with LSB disk s. The last galaxy 
in our list, F568-6 (or Malin 2), has been observed by ISprayberry et al.l (|l995h and is also bulge 
dominated. The properties of individual galaxies are summarized below. 

UGC 1378: This galaxy has a pr ominent bulge and a diffuse stellar disk. It is classified as an LSB 



galaxy with an active nucleus by lSchomber] (119981). The bulge s hows diffuse x-ray emission, pos 



sibly associated with an old stellar population (IDas et al. 



2009bl) . 



UGC 1922: The disk of this galaxy is featureless but has a bright bulge that also shows diffus e 



x-ray emission. It is classified as an LSB galaxy with an active nucleus by 



Schombert! (Il998h 



UG C 1922 is also one of the few LSB galaxies that have a significant concentration of molecular 



gas (|0'Neil & Schinnerei 



20031) 



UGC 3968: Not much is known about this LSB galaxy except that it has a prominent bulge with a 
faint disk. The 2MASS image reveals a bar assoc i ated w ith the bulge and two faint spiral arms. It 
is also classified as an LSB galaxy by lSchomberti (|l998[) but it is not clear whether the galaxy has 
an AGN. 

UGC 4219: Not much is known about this giant LSB galaxy either. According to 



Schombert 



19981) . the galaxy has a large bulge and an AGN. The LSB disk shows faint spiral arms. 
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UGC 6614: This is a well studied giant LSB galaxy (|de Blok et alJll995a). It is close to face 



on in morpholog y and has a large bulge surrounded by a ring like f eature (IRahman et al. 



Hinz et al. 



20071). The disk has faint but tig htly wound 



bulge hosts an A GN that is bright at optical dSchombert 



spira 



arms ([Pickering et al 



(INaik et al. 



19981), radio (IDas et al 



2001 



1997 ). The 



2006h and x-ray 



2010f) wayelengths. 



UGC 6754: This galaxy has an L SB disk (ISchombert & Bothunlll987r) and a prominent bulge but 
does not appear to ha ve an AGN (ISchombertil 19981) . The disk has flocculent spiral arms and only 



patchy star formation (jAmram et al 



1994). 



UGC 6968: Not much is know n about this galax y but it is described as an LSB galaxy haying a 
prominent bulge and an AGN (|Schomberti 1 1 99 8|) . There are two faint spiral arms extending out 
into the disk. 

UGC 7357: This galaxy has an LSB disk (IMac Arthur et al.i2003[) and a bright bulge. But does not 
appear to have an AGN ISchombertil 1 998b . The disk is fairly featureless. 



F568-6 (Malin 2): This is also a relatively well studied gi ant LSB galaxy. 



(Schombert & Bothun 1987), prominent bulge and an AGN (ISchombert 



t has an LSB disk 



19981). Like UGC 1922, 



it is one of the rare LSB galaxies that have a significant mass of molecular gas (|Das et al. 



2010 ). 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 HCT Data : 

The LSB galaxies were observed using the 2m Himalayan Chandra Telescope (HCT) at the Indian 
Astronomical Observatory (lAO), Hanle, which is remotely controlled from the Centre for Re- 
search and Education in Science and Technology (CREST), Indian Institute of Astrophysics (IIA), 
Bangalore. The spectra were obtained using a 11' x r.'92 slit (#1671) in combination with a grism 
#7 (blue region) and grism #8 (red region) which cover the wavelength ranges of 3700-7200 A 
and 5500-9000 A with dispersions of 1.46 A pixel"^ and 1.26 A pixel"^ respectively. The spectral 
resolution is around ~ 8.7 A (398 km s~^ FWHM or cr = 169 km s'' at Ha) for grism #7 and ~ 7 
A ( cr = 136 and 103 km s"^ at Ha and Can Triplet respectively) for grism #8. The slit was placed 
at the centre of the galaxy covering a central region of ~ 2" x 5" (1" corresponds to 415 pc at a 
redshiftof ~ 0.03). 

Data reduction was carried out using the standard tasks available within IRAfQ which includes 



Image Reduction & Analysis Facility Software distributed by National Optical Astronomy Observatories, which are operated by the Association 
of Universities for Research in Astronomy, Inc., under co-operative agreement with the National Science Foundation 
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bias subtraction, extraction of one dimensional spectra, wavelength calibration using the ferrous 

he wavele ngth calibrated spectra 



Okd (|l990h observed on the same 



argon lamp for grism #7 and ferrous neon lamp for grism #8. 
were flux calibrated using one of the spectroscopic standards of 
night and then corrected for the redshifts of the galaxies. Flux calibrated spectra were corrected for 
galactic extinction using ISchlegel et al.1 (|1998r) . The spectra were not correc ted for intrinsic dust 



extinction beca use LSB galaxies are known to have intrinsically less dust ( |Greene & Ho 



2007 



Mdetal 



20091). The blue and the red spectra were combined together with the help of scombine 



within specred package using a suitable scale factor estimated at the flat continuum portions of 
the overlapping part of the spectra. A log of the observations are presented in Table El and the flux 
calibrated HCT spectra are plotted in Figure [T] We have used the HCT observations to identify the 
LSB galaxies that host AGN activity. The black hole masses and bulge velocity dispersions were 
determined using SDSS data for all but one of these galaxies. 



3.2 SDSS DR7 Data : 

Our pilot project on LSB galaxies was started in the year 2006 wherein careful selection of objects 
for which SDSS data were unavailable was considered. But when the DR7 data was released to 
the astronomy community, we found that a few galaxies from our sample were included. The 
resolution of the SDSS data is better (~ 70 kms"'). Due to this, we have used SDSS DR7 data 
of our sample LSB galaxies for modeling/estimating the BH masses and velocity dispersions. For 
measuring the stellar velocity dispersion, SDSS ofl'ers a set of about 32 spectra of giant G and 
K stars of old open cluster M67. The stellar templates were observed in an identical manner as 
the SDSS LSB spectra and hence efl^ects arising due to template mismatch are minimal. SDSS 
spectra are observed through a fiber of 3" diameter which transforms into an area of 2.9 kpc 
diameter for a redshift of z~ 0.021. This area includes sufficient stellar light, nevertheless with 
stellar templates observed in the same setup, most of the stellar light would be removed after 
decomposition. In comparison, HCT spectra cover an area of ~ 5.4 x 13.5 kpc at the same redshift 
and the stellar templates were not observed in the same setup. Thus, for the estimation of BH 
masses after decomposing the broad and narrow components in the Ha region, and for the stellar 
velocity dispersion, we use SDSS spectra to give a conservative estimate of the above parameters. 
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4 SPECTRAL DECOMPOSITION 



The emission lines appear weak in our sample, but the Ha emission line is clearly present in several 
of the galaxies. In order to isolate the Balmer and major nebular emission lines better, we have 
decomposed the observed spectrum into different constituents. Major contribution to the observed 
spectra arises f rom the underlying s tellar population, and we use high-resolution model spectra 
of Starburst99 (jLeitherer et al.lll999h after degrading it to the resolution of the observed spectra. 
While a composite spectrum based on stellar spectral libraries would have been more realistic, 
Starburst99 is simpler to use. The decomposition was executed only up to the wavelength of 7000 
A as the Starburst99 high resolution model spectra are available only upto the above mentioned 
wavelength. 

The [Nil] A 6584 and [Oiii] A 5007 lines along with Hor a nd WIS lines are used t o calculate 



Pettini & Pagel 



(l2004h . The oxy- 



oxygen abundance assuming the empirical relation obtained by 
gen abundances and hence met a llicitie s are close to solar in value, particularly for F568-6 and 
UGC 6614 (see also UcGaugh ( 1994 ^ for UGC 6614). The exception is the galaxy UGC 7357 
which may have metallicity slightly less than s olar. We have hen ce adopted solar abundances in 



applying the Starburst99 model. The method of 



Mei et al. 



(|2009[) was used for decomposition of 



stellar light, along with emission from the Fe i an d Feu complexes, and in a few cases, a power 
law component. We use IVeron-Cetty et al.l (120011) spectra of I Zw 1 to model the Fe i and Fe ii 



complexes. Levenberg-Marquardt algorithm dPress et al. 



19931) was used for the decomposition. 



Figure |2] shows the decomposed spectra plotted at the bottom of each plot. These spectra now 
show only the gas emission due to star formation and/or the active nucleus. 

According to lSchomberti (|l998|) . in their sample of LSB galaxies, 95% showed nuclear emis- 
sion. Our sample consists of bulge dominated LSB galaxies and is a subset of the Schombert 
sample. Out of the 9 galaxies observed, we find that only 4 galaxies show broad Ha profiles emis- 
sion, along with strong emissions from [Nii], [S ii], [Oiii] and 0, which hint the presence of an 
AGN. The line fluxes of various lines are presented in Table [31 We have adopted SDSS spectra 
for estimating the fluxes of broad and narrow components of Ha and hence in the calculation of 
BH masses, as well as stellar velocity dispersions. The SDSS spectra were taken through a fiber 
aperture of 3 arcsec in diameter (corresponding to 2.9 kpc at a redshift of 0.05). The broad and 
narrow components of Ha are seperated using the fitprofs task of IRAF. Figure [3] shows the fits 
to the Ha line profiles. The broad Ha line in AGN spectra are generally asymmetric. The above 
procedure of multiple Gaussian fit results in slightly larger errors due to this. The Ha fluxes for 
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broad and narrow components, Ha luminosity, full width at half maximum (FWHM) of the broad 
Ha lines are all listed in Table [3l 

We have adopted the penalized-pixel fitting (pPxF) algorithm of ICappellari & Ernselleml (120041) 



for recovering the ste llar velocity dispersion. pPxF is gauss-hermite parametrization ( 



van der Marel & Franx 



1993 . 



Gerhard 



19931) that works in the pixel space. The reason is that in pixel space, it is easy to 



mask gas emission lines or bad pixels from the fit and the continuum matching can be directly 



taken into account (ICappellari & Emsellem 



2004|) . Also, the estimation of measurement errors are 



simplified. PPxF creates an algorithm, wherein, initial guesses for V (redshift, z) and cr are pro- 
vided. The model spectra are convolved with a broadening function using initial cr values, is 
calculated for each dataset. The residuals or;^'^ for each of the data points is perturbed and fed into 
a non-linear least squares optimization routine, in this case, monte-carlo optimization. The whole 
procedure is iterated to obtain V, cr and gauss-hermite polynomials. The estimation of gauss- 
hermite polynomials becomes a problem if the observed velocity dispersion is of order 2 pixels 
or cr < 140 km s' ^ (1 pix = 70 kms~' for SPSS d ata). A detailed explanation of the procedure can 



be obtained from 



Cappellari & EmsellemI (l2004h 



With the availability of libraries with high spectral resolution stellar and galaxy spectra, tem- 
plates can be carefully matched with the observed galaxy spectra. The 32 G and K giant stars of 
old open cluster M67 observed with the SDSS are used as stellar templates. Empirically, early K 
giants consistently provide the closest match to both the Mgib and Can Triplet regions of many 
AGN samples. Here, we have tried to measure the stellar velocity dispersion by fi tting the Can 



Greene & Ho 



Triple t lines for these galaxies UGC 6614, UGC 6968 and F568-6. According to 
(120071) . the optimal spectral region for measuring cr depends on the Eddington ratio, continuum 
level of the AGN, and redshift of interest which suggests that for z < 0.05, Ca ii Triplet is the region 
of choice for Eddington ratio < 0.5 for the most reliable measurement of cr. The pPxF code was 
applied to only Can Triplet region to estimate cr. Figure |4] displays the fits to the data; fitting error 
of ~ 10% is estimated from the residuals. 



5 AGN ACTIVITY AND BLACK HOLE MASSES 



5.1 Emission Line Diagnostic Diagram 



It is interesting to investiga te the positions occu pied by these LSB galaxi es in the BPT diagnos 



tic diagram (first given b y 



1981 



Kewlev et al. 



2001 



Baldwin et al 



Kewlev et al. 



1981 



and improved further by 



Veilleux & Osterbrock 



20061 and references therein). Figure [5] shows the diagnos- 
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tic diagrams plotted for [Oiii]/H/3 vs. [Nii]/Hq', [S ii]/HQr and 0/HQr respectively. Also plotted in 
these diagrarns are t he demarcation lines between starbursts, Se yferts and LINERS obtained from 



Kewley et al. 



200 IL 



KaufFmann et al 



2003 



Kewley et al.ll2006L From the plots, four galaxies ap- 



pear strong candidates for AGN, showing LINER-like activity. Of these, the galaxies F568-6 and 
UGC 6614 also show a high value of [Om]/HjS~1.5 when compared to other galaxies in the sam- 
ple and could be Seyfert-like. It may be not ed that the [O ml/Hjg ratio for low-mass AGNs selected 



from SDSS is 1.9 for one kind of sample (IGreene & Ho 



2007 ). The width of the broad Ha lines 



observed in these galaxies is below ~ 2000 kms ' indicating that they both belong to the class of 
narrow-line Seyfert Is (NLSls) galaxies. 



5.2 Stellar Velocity Dispersion cr* 



T he stellar velocity dispers i on cr» is measured using the Ca ii Triplet lines adopting the pPxF code 



of 



Cappellari & EmsellemI (120041 ). We could detect the Can Triplet lines at 8542 A and 8662 A 



in some of our galaxies such as UGC 6614 and F568-6. Hoever, due to template mismatch, we 
could not estimate the cr» using the HCT spectra. Hence we used the SDSS spectra for measuring 
cth.. The spectral resolution of SDSS spectra are about 4 A, which amounts to a velocity width of 
about 140 km s^ and cr = 70 km%~\ The derived values for UGC 6614, UGC 6968 and LSBC 
F568-6, are in the range 150-210 km s^^ . The SDSS spectrum is not available for the galaxy UGC 
1922 and hence we could not estimate the stellar velocity dispersion for this galaxy. We estimate 
(T, ~ 157, ~ 196 and ~ 209 kms'^ for UGC 6614, UGCC 6968 and f568-6, respectively using 
SDSS spectra. These values are shown in Table HI Figure |4] shows the fits to the observed data 
obtained using the pPxF code. 

There is some discussion in the literature on the need to reduce the observed stellar velocity 
dispersions to a uniform system, since the observed values are averaged over the slit or aperture 
size which will translate into diffe rent sizes on the face of the galaxy with respect to the bulge 
scale length. Ijorgensen et al.l (| 19951) transform the values to the equivalent of an aperture of radius 
rj%, where is the efi'ective bulge radius. The effective bulge radius is about 4'.'2 for U GC 6968 as 



calcul ated by 



Gavazzi et al. 



(|2000|) using near-IR H band image of the galaxy. Though iMcGaugb 



19941) have attempted only to fit the disk, t he bulge is visible in their surface brightness profile 



plots for the galaxies UGC 6614 and F568-6. 



de Blok et al.. (J 995b) have also obtained the surface 



brightness profile for UGC 6614 and fit only the disk. The effective bulge radius r^ for the galaxies 
UGC 1922, UGC 6614 and F568-6 are not available in the literature. On the other hand, as pointed 
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out by iFerrarese & Merrid (|2000l) . the applie d corrections for the veloc ity dispersion cr^ are very 



small, the maximum correction being < 5%. 



Greene & Hoi (I2005h and 



Pizzella et al. 



(120041) find 



that radial dependence of cr, is flat with less than 7% correction for early-type galaxies. Following 
these arguements, we have not applied any correction based on r^ to cr^,, the spectra are extracted 
from a region of 2" x 5". 



5.3 The Mass of Central Black Hole (Mbh) 



The blackhole masses are calculated using the equation given in iGreene & Hoi (l2007h using Ha 
luminosity and FWHM. The masses are ~ 0.3 x 10*^ Mq for the galaxies UGC 19 22, UGC 6968 



Greene & Ho 



and F 568-6, and lie on the lower mass tail of the low-mass blackhole sample of 
(|2007() which has a median mass of 1.3 x 10^ Mq. UGC 6614 has a slightly higher BH mass of 
3.8 X 10^ Mq. The average Lboi/LEdd ratio calculated for their sar nple is about 0.4 and suggests 



that their sources are radiating at high fraction of Eddington limit (|Greene & Ho 



2007). We find 



the values of 0.18, 0.023, 0.046 and 0.106 for UGC 1922, UGC 6614, UGC 6968 and F568- 



6 resp ectively. While these are lower compared to the median for the sample of .Greene & Ho 
(120071), they are within their observed range. 



We could identify clear signature for AGN in 4 out of 9 objects in our sample which agrees 



with the high (50%) occurrence of AGN found by 



Schombert (1998) in LBGs. 



5.4 Interesting case of UGC 6614 

The emission line spectra of UGC 6614 obtained after the decomposition of stellar light shows an 
interesting feature. A bump is noticed at the blueward of the Ha emission from the galaxy. This 
blue bump is noticeable in the observed flux calibrated spectra before spectral decomposition but 
is clearly seen after the decomposition. The blue bump is the excess emission at Ha which could 
be arising from ionised gas travelling at speeds ~ 3600 km s'^ towards us, centered at 3920 km s"^ 
The blue bump is also noticed blueward of H/3 and the velocities with which the gas streaming out 
towards us ~ 3600 km s~^ (median value of Ha and H/3) centered at 3360 km s"^ from H/3, similar 
to Ha. This emission a t Ha an d H/3 wavelengths overlap indicating the feature to be real as shown 
in Figured Das et al.l (2009a') detected a compact core and a one-sided radio jet in UGC 6614 
from 610 MHz map. A n extended feature is also indicated in a low resolution VLA map at 1420 



MHz (Das etal. 



2009a|) . The blue shifted ionized gas emission could indicate a jet or hotspot along 
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the line-of-sight. Similar asymmetric blue bump of th e [0111] lines were de tected from a bunch of 



Type 1 - Type 2 Seyferts from SDSS DR2 sample by iGreene & Hoi (|2005h . These wings indicate 
radial motions in the NLR, associated with an outflow. The outflowing components are principally 
responsible for impa rting supervirial mo tions to the gas and originate from a more compact region 



closer to the centre ioreene & Ho 



20051) 



5.5 The M 



BH 



(T, Plot 



The Can Triplet line widths and masses of blackholes for the three galaxies in our sample are 
shown in the Mrh - cr^ plot in Figure [71 Also plot t ed in the figure are the line a r regr ession 



lines given by 



Giiltekin et al 



(l2009h . 



Tremaine et al. 



(l2002h and 



Ferrarese & FordI (l2005h (dot 



ted, dashed and solid lines respectively) for Mbh against cr^. The low-mass AGNs hosted within 
LSB galaxies occupy the region ju st below the lowest mass blackhole of Circinus galaxy from 



the sample of 



Giiltekin et al. 



(|2009|) . well belo w extrapo 



other hand, three AGNs in LSBs observed by 



2.8 - 20 X 10^ M„ lie closer to the 



Mei et a 



Tremaine et al 



ations of high-mass blackholes. On the 
(|2009h in the blackhole mass range of 



(|2002h relation, thou gh systematica ' 



would be of interest to study more LSBGs and low luminosity AGN of 
better understanding of faint luminosity end of Mbh - cr* relation. 



Greene & Hoi ( 2007 ) for a 



ly low er. It 



6 DISCUSSION 



1. Intermediate Mass Black Hole (IMBH) in GLSB galaxies: One of the main results of our 
spectroscopic study is the detection of broad Ha emission in GLSBs and the subsequent estimation 
of nuclear black hole masses from SDSS spectra in bulge-dominated GLSB galaxies. The AGNs 
fall in Seyfert- LINER region in the diagnostic diagram (refer figure [5]). 

We obtained masses ~ 3 x 10^ Mq for three GLSB galaxies in our sample, which fall in the 
IMBH range rather than the SMBH range. A higher blackhole mas s is estimated for UGC 6614 



which is ~ 3.8 x 10^ M© and a simil ar estimate was given ea rlier by 



Das et al 



low resolution optical spectrum from lSprayberry et al. 



(I2009bt) based on a 
(Il995h . Another estimate of ~ 1 .2 x 10^ Mq 



Naiketal 



(|2010|) . which is lower than 



was derived later from AGN X-ray variability studies by 
the present estimate. It must be borne in mind that these estimates are fairly approximate as they 
are based on the as sumption that the gas in the broad line region in the AGN is in virial equilibrium 



Kaspi et al 



2000 ). 



IMBHs are fairly rare in the galactic nuclei and have been detected mainly in late type spi- 



rals (iFilippenb 
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o & nlioos 



Greene & Ho 



2004 



Satyapal et al.ll2007l) . nearby galaxies (|Seth et al. 



20101) or dwarf galaxies (|Barth et alJ 120041 ). They are difficult to detect dynamically at large dis- 
tances; hence AGN activity is one of the main methods through which we detect them. The pres- 
ence of IMBHs in GLSB galaxies is surprising as their bulges are well developed; in fact a SMBH 
would be far more typical for these bulges. This suggests that the lack of disk evolution in these 
extreme late type galaxies has affected the evolution of their nuclear BHs. Galactic disk activity 
contributes to the growth of SMBHs through gas inflow, star formation and mass accumulation in 
the nuclei of spiral galaxies as observed in bulge dominated, star forming early type spirals. Large 
scale disk instabilities such as bars and spiral arms exert gravitational torques that funnel gas into 
gala xy centers leading to n uclear star formation and the build-up of central mass concentrations 



(e.g. 



Friedli & Benz 



1993h. This can r esult in the growth of nuclear black holes and bulges in 



galaxies (I Kormendy & Kennicutti 120041) . This process of disk evolution leading to the growth of 
central mass concentrations i n galaxies is prevented fro m happening or slowed down when there is 



a dominant dark matter halo ( Ostriker & Peebles 



1973 ). Thus, the lack of disk evolution and rela- 



tively low mass of the black hole may share the same origin - which is the presence of a dominant 
dark halo in the galaxy. 

2. Constraining the M - cr relation for extreme late type galaxies: In the past ten years 
the M - cr and M - L relations have become established benchmarks for galaxy and BH evolu- 
tion theories. Our present work affects these correlations in two ways; first it helps constrain the 
ow mass end of the M - a correlation and secondly it helps to constrain the scatter in the plot 



( Gultekin et al. 



20101) . The low mass end of the M - cr relation is populated by late type galaxies 



or dwarf galaxies; many are often outliers in the plot. Our present work shows that extreme late 
type galaxies are also fairly offset from the main M - cr line (Figure |7]). It also suggests that dwarf 
galaxies and extreme late type galaxies have different evolutionary paths compared to early type 
galaxies and the more massive ellipticals at the high SMBH end of the M - cr relation. Models of 
galaxy evolution thus need to incorporate late type systems such as GLSB galaxies in their overall 
picture. 



The late type spirals also increase the scat ter in the corre 



ellipticals and early type spirals are included (|Tremaine et al. 



ation which is tighter wh en only 



2002 . 



Beifiori et al. 



20091) . Many 



theoretical stu dies have been undertaken to expl ain the M-cr correla tion and predict the high mass 



end of the plot (INatarajan & Treister, 20091) . iDalla Bonta et al.l (|2009l) carried out HST observations 
of three brightest cluster galaxies (BCGs) and estimated masses of SMBHs to be ~ 10^ Mq present 
in these BCGs. While for one galaxy, SMBH mass correlates well in the M-cr and M-L plot at the 
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high mass end, the other two galaxies show inconsistencies with the two relations. The sample is 
small to derive any conclusio ns, but hints that there co uld be scatter in the SMBH scaling relations 



at the high mass end as well (|Dalla Bonta et al. 



2009 ). 



However the low rtiass e nd does not appear to have a clear cutoff according to most models 



(|Volonteri & Natarajanll2009l) . In fact, we could be missing observationally a large fraction of the 



lower mass BHs in the centres of galaxies. Thus there is an increase in the scatter at the low mass 
end of the correlation. This scatter could be due to m easurement erro rs or could be intrinsic to 
the BH evolution processes in the galaxies, themselves (I Volonterill2007|) . Our present study is thus 
important for understanding the overall trends in the low mass end of the BH mass spectrum. 

3. AGN evolution in late type galaxies: Studies have shown that the space density of high 
luminosity AGNs peak at redshif t of z ~ 2; this is also the redshift at w hich the most massive 



SMBHs were formed in galaxies (|Cowie et al 



2003 



Hasinger et al 



2005 ). However, in the local 



universe the most rapi dly growing BHs appear to be those in the lower mass range of 10^ - 10^ 



(|Goulding et al 



2010h . Also, studies of nearby galaxi es show that it is the mo st massive BHs in late 



type galaxies that are growing at the present epoch (|Schawinski et al.ll2010h . The GLSB galaxies 



in our sample fall into the latter category as they have large bulges; though they have lower BH 
masses, they appear to be accreting and hence luminous in the optical domain. 



Das et al. 



(l2009bh . 



4. Decoupled Bulge-Disk Evolution in GLSB galaxies: As suggested by 
the bulges of GLSB galaxies appear to be very evolved compared to their disks. In general bulges 
form in two ways; one is through repeated galaxy mergers or accretion events that lead to the 
formation of a central spheroidal mass distribution dSpringel et al.ll2005h . The second is through 
secular evolutionary processes where disk instabilities lead to bars, spiral arms, gas infall and 
the evolution of a disky pseudobulge (|Kormendy & Kennicutti 120041) . These processes result in 
disk star formation and enhanced disk structure, both of which are not observed in most GLSB 
galaxies; instead their disks are metal poor and often fairly featureless. So the bulges in GLSB 
galaxies probably formed in a different way; one possibility is th at galaxy mergers resulted i n 



spheroidal bulges and then the disks were rebuilt from accreted gas (ISpringel & Hernquist 



2005 ). 



Such an evolutionary scenerio would lead to a bulge that is relatively decoupled from its disk or 
its central black hole. 

5. BHs in Halo Dominated Galaxies: Although the correlation of BH mass and galaxy prop- 
erties is now we ll estab lished, it is still not clear exactly what regulates black hole growth (e.g. 



Booth & Schayd (|201ol) and references therein). Mass accretion close to the black hole, bulge 



mass and the mass of the dark matter halo are some of the factors important for regulating black 
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hole growth in galaxies. It is not clear which factor is the most important or whether all the 
processes play a role. Several theoretical studies have e xplored how the potential of the dark 



halo may regulate bulge evolution and black hole growth (|Booth & Schaye 



Silk & Rees 



2010 



Xu et al 



200' 



19981) and there are observationa 



halo and black hole mass (IBaes et al. 



2003 



studies that indicate a correlation between the dark 



Ferrarese 



20021 : 



Pizzella et al, 



20051) . On the other 



hand, some studies show that nuclear black holes masses do not correlate with the dark halo mat- 
ter haloes of galaxies and dark matte r gravity is not directly responsible for black hole growth 



Kormendy & Bender 



2011 



Ho 



2007h . Such studies suggest that SMBHs co-evolve with classical 



bulges or ellipticals only. GLSB galaxies are halo dominated and often bulge dominated as well. 
Hence they are ideal systems to study the dark halo-BH relation and this should be investigated in 
future studies. 



7 CONCLUSION 

(i) The paper presents spectroscopic observations of the nuclear regions of nine low surface 
brightness galaxies observed in the wavelength range 3700-9000 A. The stellar light has been 
subtracted from the nuclear spectra to obtain only gas emission spectra. Broad Ha lines along 
with strong [Nii], [Sii], [Oi] lines are detected in four galaxies namely UGC 1922, UGC 6614, 
UGC 6968 and F568-6 confirming the presence of AGN activity in these LSB galaxies. SDSS 
spectra of three galaxies namely, UGC 6614, UGC 6968 and F568-6 are used to estimate the BH 
masses and stellar velocity dispersion as the resolution of the SDSS spectra is good enough (~ 70 
km s~') to give a conservative estimate on the above parameters. 

(ii) The BPT AGN diagnostic diagram was created using the emission line ratios. It is clearly 
seen that the above four galaxies lie in the AGN regime and more closely, in the Seyfert regime. 
UGC 3968 might also host a starburst-AGN composite nucleus. 

(iii) The broad Ha line widths (900-2500 kms~^) and luminosities (10^^ erg s"') are used to 
deduce the nuclear blackhole masses in the aforementioned galaxies; the masses for three galaxies 
are ~ 3 x 10^ Mq and for UGC 6614, the BH mass is estimated to be about 3.8 x 10^ Mq. The 
masses suggest that the nuclei of LSB galaxies have IMBHs rather than the SMBHs found in the 
centres of brighter galaxies. UGC 6614 also shows an interesting feature of a blueshifted bump of 
He emission which can be attributed to outflow of gas travelling at speeds of 3600 km s"^ towards 
us. The blue bump feature is detected in the H/3 region as well. There could be more such AGN 
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in the sample that may be identified through other means such as the reverberation technique with 
improved sensitivity. 

(iv) The stellar velocity dispersion, cr* is measured for the three galaxies UGC 6614, UGC 
6968 and F568-6; the values lie between 150-210 km s"^ The three low-mass BHs lie below the 
standard line in the M - cr» plot and lower than the ones studied by 



Mei et al 



mm 
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Figure 1. Spectra of the LSBs observed from HCT using the blue (#7) and red (#8) grism and combined together using the scombine 
task of IRAF. The x-axis represents the rest-frame wavelength. The spectra are arbitrarily displaced for clear viewing. 



Table 1. The sample of LSB galaxies 



Galaxy 


Other 


Galaxy 


Vel." 


Distance* 


Galaxy 


Galaxy 


Galaxy 


Absolute 




Name 


Type 


kms"' 


(Mpc) 


Coordinates (J2000) 


inclination 


size (arcmin) 


Magnitude (Mh) 


UGC 1378 


PGC 007247 


(R)SB(rs)a 


2935 


43 


01h56ml9.2s +73dl6m58s 


68.6« 


3.4 


-24.01 


UGC 1922^ 


PGC 009373 


S 


10894 


151 


02h27m45.8s +28dl2m33s 


29.6" 


2.1 


-25.25 


UGC 3968 


PGC 021636 


SB(r)c 


6780 


94 


07h42m45.2s +66dl5m30s 


37.8" 


1.4 


-23.55 


UGC 4219 


PGC 022766 


SA(rs)b 


12433 


170 


08h{)6m42.8s +39d05m25s 


44.7" 


1.1 


-24.37 


UGC 6614 


PGC 036122 


(R)SA(r)a 


6352 


86 


Ilh39ml4.9s +17d08m37s 


29.9" 


1.7 


-23.70 


UGC 6754'' 


PGC 036740 


SA(rs)b 


7025 


96 


Ilh46m47.2s +20d40m32s 


29.1" 


3.0 


-24.74 


UGC 6968 


PGC 037704 


S 


8232 


113 


Ilh58m44.6s +28dl7m22s 


44.7" 


2.8 


-24.04 


UGC 7357 


PGC 039640 


SAB(s)c 


6682 


91 


12hl9ml3.4s +22d25m54s 


48" 


1.6 


-22.17 


F568-6 


Mahn 2 


Sd/ 


13830 


189 


10h39m52.5s +20d50m49s 


38" 


1.5 


-25.50 



" Heliocentric velocities from NED Galactocentric distances from NED assuming Ho = 73 kms" 'Mpc ^ = IC 226 = NGC 3883 



Xu, B., Wu, X., & Zhao, H. 2007, ApJ, 664, 198 
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Figure 1. - Continued 



Table 2. Observation details of long slit spectra obtained from HCT. 



Galaxy 


Date of Obs. 


Exptime" 






Gr7 / GrS in sec 


UGC 1378 


22-11-2006 


3600 / 3600 


UGC 1922 


22-12-2006 


3600/3600 


UGC 3968 


22-11-2006 


1800/3600 


UGC 4219 


19-02-2007 


2400 / 2400 


UGC 6614 


05-07-2006 


3600/3600 


UGC 6754 


19-02-2007 


1800 / 2400 


UGC 6968 


14-05-2007 


1800/2400 


UGC 7357 


19-02-2007 


2400/1510 


F568-6 


22-12-2006 


3600/3600 



" - Exposure times in the grism 7 (blue) and grism 8 (red). 
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Figure 2. Plots show observed HCT spectra on the top, wit h best fitting Gyr model sh owing the age of the underlying stellar population 
as a dashed line. Fe i and ii template spectra taken from IVeron- Cettv et alj bOOlh are plotted as dot-dashed. Dotted lines seen in 
UGC 1378, UGC 3968 and in UGC 4219 are power-law continuum following lMei et alj i2009l) . The bottom spectra drawn as thick 
solid lines show the subtracted spectrum. The x-axis represents the rest-frame wavelength. 
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Figure 2. - Continued. UGC 6614 and UGC 6968 are best fit using only the stellar spectra from Starburst99, while the spectra of UGC 
7357 and F568-6 require additional power-law continuum shown as dotted lines in this figure. 
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Figure 3. Fits to the broad and narrow emission line components for the galaxies F568-6 and UGC 1922 from SDSS. The points show 
the observed spectra after decomposition. The dotted lines show the individual gaussian fits to the lines and the solid line shows the 
combined fits for the region. The spectra used for estimating BH mass in F568-6 is the SDSS spectra while for UGC 1922, we have 
used HCT spectra as SDSS data for this galaxy is unavailable. 
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Figure 3. - Continued using the SDSS spectra for the two galaxies UGC 6614 and UGC 6968. 
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1- UGC 1922; 2- UGC 6614; 3- UGC 6968; 4- LSBC F568-6; 5- UGC 7357; 6- UGC 3968; 7- UGC 4219. 



Figure 5. Diagnostic diagrams for the sample of LSB galaxies, observed f rom HCT. The dashe d line in the leftmost figure, 
log([0 in]/Hj8) vs log([N iil/Hff), is the starburst-Seyfert demarcation taken from Kaulf mann et al The dotted line is the extreme 

starburst line of Kewley et al. 1 2001). The solid lines showing the demarcation in log([0 iii]/H/3) vs log([S ii]/Hq;) and log([0 in]/Hy3) vs 
log([0 il/Ha) i s taken fromlKewlev et al.N200l]) . The straight line drawn in the centre and right plots is the Seyfert-LINER demarcation 
line taken from lKewlev et al.H200d) . The [Oi] emission was clearly detected only in four of the galaxies UGC 1922, UGC 6614, UGC 
6968 and F568-6. These galaxies definitely host near-Seyfert kind of nuclei. UGC 7357 and UGC 4219 clearly fall in the starburst 
regime, and the galaxy UGC 3968 might be hosting a starburst-AGN composite nucleus as seen from the diagnostic diagrams. 
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Figure 6. Line fits for the region around Ho- and H/3 for the galaxy UGC 6614. The x-axis is interms of velocities. A blue bump is 
clearly noticeable around Ha and HyS which is centered at -3600 ± 300 km s"' and moving at a velocity of ~ 3600 km s"' towards us. 




Figure 7. M-cr plot for the LSBs. Points represe nt the Mbh and velocity d ispersion of the 3 LSBs UGC 6968, UGC 6614 a nd F568-6. 
The solid line is t he Mb/z-ct, relation ta ken from lFerrarese & Merritll feOOOl) . the dashed line is from lTremaine et al] j2002h and dotted 
line is taken from lGultekin et alj <2009l) . The Mbh and cr, are estimated from the SDSS spectra. 



to 

Table 3. Fluxes of emission lines in the units of 10 " erg cm - s ' obtained from HCT spectra. <^ 



Galaxy 


[Oil] 




[OIII] 


[OI] 


Ua" 


[Nil] 


[SII] 




3727 A 


4861 A 


5007 A 


6300 A 


6563narraiv 


6548 A 


6584 A 


6717 A 


6731 A 


UGC 1378 




** 


5.75 ± 0.54 




5.30 ± 0.81 




8.20 ± 0.87 


5.30 ± 1.56 


4.74 ± 1.25 








(9.41 ± 0.84) 




(10.12 ± 1.19) 




(10.16 ±0.98) 


(13.07 ±3.12) 


(12.72 ±3.03)) 


UGC 1922 


31.01 ± 0.07 


4.35 ± 0.67 


3.32 ± 0.56 


2.21 ± 0.73 




3.77 ± 0.69 


8.29 ± 0.59 


5.20 ± 1.36 


3.56 ±0.89 




(9.60 ± 0.23) 


(9.47 ± 1.08) 


(8.44 ± 1.23) 


(9.82 ± 3.06) 




(15.58 ± 1.88) 


(18.72 ±3.08) 


(11.67 ±2.60) 


(9.93 ± 2.95) 


UGC 3968 




2.08 ± 0.13 






3.50 ± 0.51 


1.81 ±0.52 


2.75 ± 0.35 


0.91 ± 0.40 


1.45 ±0.59 






(8.76 ±0.71) 




** 


(10.41 ± 0.92) 


(14.73 ±3.01) 


(9.47 ±0.91) 


(8.02 ±2.75) 


(13.31 ±4.58) 


UGC 4219 




1.10 ± 0.25 




** 


4.45 ± 0.50 


0.58 ± 0.35 


1.08 ±0.39 


0.69 ± 0.45 


0.46 ± 0.38 






(10.03 ±2.59) 






(11.75 ± 0.86) 


(8.91 ±5.37) 


(10.30 ±3.76) 


(11.37 ±9.13) 


(10.64 ± 10.57) 


UGC 6614 


7.89 ± 0.09 


2.56 ±0.35 


5.12 ±0.31 


3.78 ±0.86 




1.69 ±0.55 


5.58 ±0.40 


5.46 ± 0.36" 


d 




(10.52 ± 1.80) 


(12.41 ±0.94) 


(13.93 ±2.88) 


(15.63 ±0.98) 


(16.78 ±3.91) 


(15.63 ±0.98) 


(18.55 ±2.13) 


UGC 6754 






0.59 ± 0.38 




1.41 ±0.83 




0.18 ±0.42 


0.54 ± 0.33 


0.76 ± 0.56 


d 






(8.59 ± 6.63) 




(5.54 ± 8.57) 




(5.83 ±4.84) 


(7.34 ± 4.46) 


(11.15 ± 8.04) 


UGC 6968 




4.52 ± 0.62 


3.49 ± 0.43 


1.14 ±0.41 




0.40 ± 0.53 


2.10 ±0.53 


1.49 ±0.61 


2.14 ± 1.09 


d 




(15.80 ±2.58) 


(12.25 ± 1.83) 


(7.45 ± 3.57) 








(6.78 ± 2.60) 


(15.06 ±7.56) 


UGC 7357 




2.03 ± 0.08 


0.45 ± 0.08 




4.60 ± 0.06 


1.16 ±0.09 


2.57 ± 0.07 


2.31 ±0.16 


1.83 ±0.13 


d 




(8.51 ±0.25) 


(3.95 ± 0.48) 




(8.21 ±0.16) 


(10.37 ±0.78) 


(9.85 ± 0.36) 


(10.25 ±0.36) 


(11.12 ±0.79) 


F568-6 




1.66 ±0.26 


1.55 ±0.31 


0.49 ± 0.39 




** 


2.34 ± 1.99 


0.79 ±0.61 


0.54 ± 0.39 


d 




(12.57 ±2.23) 


(10.17 ± 1.56) 


(9.65 ± 2.33) 




** 


(15.31 ± 1.31) 


(11.56 ±7.78) 


(9.77 ±7.51) 



9^ 



" - Hff fluxes for the objects which do not host AGN. 

- broad Ha line is detected and hence separation of broad and naiTow lines is done. These are deblended and given in the Tableg] 

- The telluric line A6S10 A exactly coincided with the [S ii] emission. After the telluric line correction, the lines could not be resolved and hence the lines could not be deblended. Hence total [S ii] flux is given here. 

- The values quoted inside the brackets, in the second row for each of the galaxies are FWHM of each of the lines whose fluxes are given above. 
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Table 4. Ha flux for narrow and broad components in erg cm"^ s"', luminosity of broad Ha Une, I^WHM of the broad Ha Une, cr, 
calculated using Can Triplet lines from SDSS spectra for three galaxies UGC 6614, UGC 6968 and F568-6. Due to unavailability of 
SDSS data for UGC 1922, HCT spectra is used for estunation of BH mass. 



Galaxy 


Ha" 

Broad Narrow 


xlC^^ergs-' 


fwhrn* 

km s~' 


km s~' 


Mbh 
xlO*Mo 


UGC 1922 


8.13 ±1.50 


4.79 ± 1.93 


2.21 ± 0.410 


855.7 ± 140.7 


** 


39+° '* 

"■■^ -0.15 


UGC 6614 


32.77 ± 3.62 


5.11 ±0.71 


2.91 ±0.020 


2456.5 ± 167.3 


157.3 ±13.1 


3 89+1-21 


UGC 6968 


3.54 ± 1.68 


1.83 ±0.69 


0.54 ± 0.001 


1244.7 ± 742.8 


195.8 ± 16.4 


45+0" 

-0.38 


F568-6 


2.01 ± 0.84 


0.93 ± 0.40 


0.85 ± 0.020 


899.6 ± 372.4 


209.1 ± 18.0 


29+0-32 



Ha fluxes of the broad and narrow components in vmits of 10 erg cm ^ s ^ 
*- FWHM of the broad line of Ha. 



